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Preamble. This series made up of two articles is devoted to a new method for the spatially differentiated assessment of the impacts of 
primary airborne pollutants on human health within Life Cycle Assessments. The first part describes the method and provides exem¬ 
plary results for site-dependent exposure efficiencies. The second part deals with the application of the method within a Life Cycle 
Assessment of natural gas vehicles. 


Abstract. A new method for the spatially differentiated assess¬ 
ment of impacts of airborne pollutants on human health is pre¬ 
sented. It is'applicable to primary pollutants with linear expo¬ 
sure response functions. This includes the most important 
primary air pollutants from transportation and energy genera¬ 
tion. The article looks at the spatial differentiation of impacts 
due to emission height and the local population density distri¬ 
bution around the emission site, as has been predicted using a 
Gaussian plume model. The differentiation due to population 
density is captured by way of five generic spatial classes: large 
cities in agglomerations, highly densified districts in agglo¬ 
merations, cities in urbanized regions, country average districts, 
and low density rural districts in rural regions. Average impacts 
are calculated for each class. The method is simple enough to be 
applied to a large number of emissions within Life Cycle As¬ 
sessments. It was used to calculate site-dependent exposure 
efficiencies for a variety of primary pollutants emitted at differ¬ 
ent heights. For traffic emissions of pollutants with short at¬ 
mospheric residence times, the exposure efficiencies vary by a 
factor of 5 across Germany and by a factor of 75 across Europe. 
This differentiation due to population density decreases signifi¬ 
cantly with an increasing atmospheric residence time of the pol¬ 
lutants and with an increasing emission height. 

Keywords: Airborne pollutants; emission height; exposure effi¬ 
ciency; human toxicity; Life Cycle Impact Assessment (LCIA); 
population density; spatial differentiation 


Introduction 

Methods for the spatial differentiation of the assessment of 
local or regional impacts within Life Cycle Assessment (LCA) 
are a topic of current interest. This applies to human health 
impacts in particular. The present paper addresses this issue 
for the case of carcinogenic and respiratory health impacts 
of primary airborne pollutants. Only atmospheric exposures 
are taken into account. For the most important air pollut¬ 
ants from transportation and energy generation, a linear 
exposure-response function can be assumed (Spadaro and 


Rabl 1999). In this case, the marginal impact of an emission 
with a given mass depends upon the emission height, as well 
as on the population density distribution and the meteoro¬ 
logical conditions (but not on the background concentra¬ 
tion of the pollutant) around the emission source. In princi¬ 
ple, it is desirable to take these spatially variable parameters 
into account in the impact assessment. However, due to the 
large number of processes to be considered in an LCA, there 
is a trade-off between the spatial detail and the ease of use 
of impact assessment methods. 

Existing methods of impact assessment are located at differ¬ 
ent positions along this trade-off. The methods commonly 
used at present are generic, i.e. they do not include any spa¬ 
tial differentiation. Examples are the method of Critical 
Volumes (Habersatter 1990), the CML method (Heijungs et 
al. 1992) and the Critical Surface Time method (Jolliet and 
Crettaz 1997). Hofstetter (1998) considers variations of the 
population density at a sub-continental scale (i.e. between 
countries or sub-continental regions) in the context of mul¬ 
timedia models, but does not include differentiations with 
regard to emission height or settlement structures at a local 
scale, i.e. between cities and rural areas. 

At the other end of the trade-off spectrum, the use of site- 
specific impact pathway analyses in the context of LCA has 
been proposed by Krewitt et al. (1998). Exemplary damage 
factors for a number of European sites are provided by 
Spadaro and Rabl (1999). A software tool is available to 
facilitate such site-specific assessments (IER 1998). How¬ 
ever, site-specific meteorological data (including local wind 
directions) still need to be added by the user. Their collec¬ 
tion might not always be feasible for a large number of 
sources or for traffic sources which pass by many different 
sites. Potting (2000) provides a framework for the calcula¬ 
tion of site-specific impacts, which is based on (IER 1998). 
Typical meteorological data for four zones within Europe 
are used, but the issue of local wind direction is not ad¬ 
dressed, and no operational guidance for the determination 
of local population densities is provided. Moriguchi and 


Int J LCA 6 (5) 257 - 264 (2001) 

© ecomed publishers, D-86899 Landsberg, Germany and Ft. Worth/TX, USA • Tokyo, Japan • Mumbai, India • Seoul, Korea 








Generic Spatial Classes for Human Health Impacts, Part I 


LCA Methodology 


Terazono (2000) also avoid the issue of local wind direc¬ 
tions by hypothetically assuming a uniform distribution with¬ 
out justification. 

The method presented here provides both a detailed justifi¬ 
cation for avoiding the tedious collection of site-specific wind 
directions and operational data for local population density 
distributions. It follows an intermediate way between a ge¬ 
neric and a site-specific approach: The differentiation of the 
health impacts between cities and rural areas is approxi¬ 
mately captured by a limited number of generic spatial 
classes, as suggested by Potting and Hauschild (1997) and 
recommended by Udo de Haes et al. (1999) (section 1). On 
the basis of generic spatial classes defined for Germany, site- 
dependent exposure efficiencies were calculated for a vari¬ 
ety of airborne pollutants (section 2). 

1 Methodology 

Consider the emission of a mass M of a primary air pollut¬ 
ant with a linear exposure response function of slope E (ef¬ 
fect factor) from a point source located at site i. Its incre¬ 
mental health impact AD; aggregated across the entire 
population can be written as 

AD; = EMAI, = EMEEyV (1) 

where AI ; is the incremental population exposure per emit¬ 
ted mass (conveniently expressed in units of persons pg/m 3 
year/kg). The product of Al ( with the volume V of inhaled 
air per time and person yields the dimensionless exposure 
efficiency EE,, which represents the fraction of the emitted 
pollutant mass inhaled by the population (Crettaz 2000). 
The exposure efficiency can be calculated by assuming a 
stationary situation as 

EE| = V/Q j r 2 J AC|(r,<p) P;(r,<p) dtp dr (2) 

o o 

where (r,cp) are two-dimensional polar coordinates around 
the site i, Ac ( (r,<p) is the stationary pollutant concentration 
corresponding to the arbitrary time-independent emission 
rate Q (note that Ac;(r,<p)/Q is independent of Q), and p ; (r,<p) 
is the population density distribution (Nigge 2000) 1 . The 
upper integration limit R needs to be determined in such a 
way that the major part of the exposure is captured. For 
pollutants with atmospheric residence times of up to several 
days, R needs to be on the order of several thousand km 
(European Commission 1995). 

In the following, it will be suitable to split the radial integra¬ 
tion interval into a short range [0,100 km] and a long range 
[100 km, R], The long-range contribution EE ; far only varies 


1 The assumption of a stationary situation can be justified as follows: 'linear ex¬ 
posure-response function' means (a) that the health impact of the pollutant is a 
function of the exposure, i.e. the time-integral of the concentration (i.e. Haber's 
law applies) and (b) that this function is linear. Since the partial differential 
equations describing the dispersion of primary pollutants are linear in both 
emission rates and concentrations (Seinfeld and Pandis 1998), the arguments 
set forth by Heijungs (1995) show that the health impact is independent of the 
time pattern of the emission (for a given set of atmospheric conditions) and can 
be calculated for the stationary case. 


on a sub-continental scale. It can therefore be approximated 
by a constant value for all emission sites within one country 
(Nigge 2000). The spatial differentiation between emissions 
in large cities and in rural areas is associated with the short- 
range contribution EE; near , which will be considered in more 
detail in the following. 

For a given emission height, Aq(r,(p) depends on the mete¬ 
orological conditions, and in particular on the frequency 
distribution of wind directions at the emission site. These 
wind roses vary on a scale of about 10 km due to local topo¬ 
graphical conditions (Zenger 1998). Site-specific meteoro¬ 
logical data are therefore very difficult to obtain. This prob¬ 
lem can be circumvented by considering classes of emission 
sites rather than the individual sites. Average exposure 
efficiencies <EE; ncar > for these classes can be calculated with¬ 
out first knowing EE, ncar for each individual site i. This is 
achieved by using a statistical argument which is illustrated 
in Box 1. 

Box 1: Statistical argument for the calculation of <EE, near > 

Consider the class of highly densified districts within agglomerated re¬ 
gions (class (1,2) in Fig. 2) as an example. This class consists of 780 
municipalities which are situated very close to the 56 largest German 
cities (class (1,1)). For an emission within one of its individual munici¬ 
palities, the wind direction therefore has a significant impact on its health 
impact: if the neighboring large city is downwind from the emission 
source (in terms of the predominant wind direction), the impact is higher 
than if the large city is upwind. Considering all municipalities within the 
class, the large city can be expected to be downwind in some cases 
and upwind in an equal number of other cases. In the class average 
impact, the effects of the wind directions at the individual municipalities 
therefore cancel each other out. 


For its mathematical formulation and empirical verification, 
see (Nigge 2000) 2 . Applying the statistical argument leads to 


<EE;> = V/Q j <Ac;(r)> <p;(r)> 2 Jtr dr + <EE i-far > (3) 
o 

The class average exposure efficiency <EE,> therefore only 
depends on the angular average Ac;(r) of the pollutant con¬ 
centration, which is independent of the wind rose at the 
emission site (Nigge 2000). This eliminates much of the site- 
specificity. Ac;(r) can be calculated on the basis of generic 
meteorological data, which may vary between countries or 
sub-continental regions, but are not site-specific. For Ger¬ 
many, such generic meteorological data (namely combined 
frequency distributions of atmospheric stability and short¬ 
term wind speeds averaged over all seasons and all times 
during the day) were generated on the basis of results from 
long-term meteorological observations (Manier 1971; Chris- 
toffer and Ulbricht-Eissing 1989). The only input variable 
required for that purpose is the annual mean wind speed. 
Three geographical zones with annual mean wind speeds in 
the ranges of 2-3 m/s, 3-4 m/s and 4-5 m/s can be distin¬ 
guished if extreme situations on the coast or in river valleys 
are disregarded (Gerth and Christoffer 1994). 

2 The statistical argument requires some modification for coastal areas. For 
Germany, which is considered here, coastal areas were disregarded since they 
only represent a small part of the country. 
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Fig. 1: Concentration ot particles (PM 10) emitted at a rate of 1 kg/h at 
different heights 


A Gaussian plume model (Janicke 1998) was used to calcu¬ 
late <ACj(r)> from the generic meteorological data for 11 
different effective emission heights h i . The model fulfils the 
requirements of the respective guideline of the German As¬ 
sociation of Engineers (VDI 1992). Fig. 1 shows radial con¬ 
centration profiles <Ac ; (r)> for particles (PM 10) emitted at 
different effective heights for the average annual mean wind 
speed of 3.5 m/s in Germany. 

Equation (3) suggests the use of the radial population den¬ 
sity distribution p ; (r) within a circle of radius 100 km around 
the emission source as the criterion for the definition of ge- 


3 The effective emission height is the sum of the stack height and the final value 
of the plume rise due to mechanical impulse or thermal effects (VD11985). In 
the following, the term emission height is used in the sense of the effective 
emission height. 
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Fig. 2: Settlement structure classes in Germany (BBR 1998a). Scale 1:7,000,000 
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Fig. 3: Radial population density distributions for the settlement structure 
classes in Germany 

neric spatial classes. For Germany, this criterion was opera¬ 
tionalized on the basis of an existing official classification 
of settlement structures (BBR 1998a). The nine settlement 
structure classes range from large cities within agglomer¬ 
ated regions to lower density rural districts within rural re¬ 
gions (Fig. 2). Each of the classes contains between 50 and 
3000 municipalities. For each municipality within a class, 
the radial population density p ; (r) was calculated in incre¬ 
ments of Ar = 10 km from data provided in (BBR 1998b) 
and then averaged across the class to yield <p ; (r)> (Fig. 3). 

The main difference between the classes is the value of <p ; (r)> 
within the first 10 km (referred to as p 10km ), w ith a spread by 
about a factor of 20 between the highest and lowest value. In 
the agglomerated regions, <p s (r)> also differs significantly be¬ 
tween classes for r > 10 km, but approaches the country aver¬ 
age population density (p D = 230 persons/km 2 ) for r —> 100 km. 
In the urbanized regions, <p ; (r)> is close to p D throughout 
the range of 10 - 100 km for all classes. Increases and de¬ 
creases of <p;(r)> in this range are due to statistical fluctua¬ 
tions. For rural districts in rural regions, p D is approached 
from below, since they are far away from more densely popu¬ 
lated areas. The fact that <p;(100 km)> = p D for all settle¬ 
ment structure classes is the reason why r = 100 km was 
chosen as the border between the short range and the long 


range, and why the long-range contribution is approximately 
the same for all sites within the country. 

By combining these radial population density distributions with 
radial concentration profiles <Ac l (r)> such as the ones shown 
in Fig. 1, the short-range contribution <EE, near > can be calcu¬ 
lated according to equation (3). The long-range contribution 
<EE| h p> was calculated with the EcoSense software (IER1998) 
as the mean value for 41 emission sites spread evenly across 
Germany. Incremental pollutant concentrations across West¬ 
ern, Central and parts of Eastern Europe are calculated with 
a Lagrangian air pollution model (Wind rose Trajectory 
Model) with a spatial resolution of 100 km x 100 km and 
combined with a database of population densities. Com¬ 
plete vertical mixing of the pollutants immediately after their 
emission is assumed, i.e. effects of emission height are ne¬ 
glected. The low horizontal resolution and the lack of verti¬ 
cal resolution of the pollutant concentration are sufficient 
for the long-range, but not for the short-range contribution 
to the exposure efficiency, which is sensitive to the emission 
height (see Fig. 1) and the variation of the population den¬ 
sity on a scale of 10 km (see Fig. 3). 

So far, only point sources were considered. This does not 
represent a limitation, however. Consider a traffic emission, 
for example, which can be represented as a line source. If all 
infinitesimal segments on the line belong to one generic spa¬ 
tial class, they are all characterized by the same <P](r)> and 
the same <ACj(r)>. Each segment therefore contributes the 
same amount to the exposure efficiency <EE,>. The fact that 
different individuals may be affected by the emissions from 
different segments does not matter because only the popula¬ 
tion exposure (rather than individual exposure) is consid¬ 
ered in equation (1). The line source can therefore be con¬ 
densed into a point source which emits the same pollutant 
mass Q per unit time. Similarly, a line source extending across 
n different spatial classes can be condensed into n point 
sources, one for each class. In that case, the share of the 
pollutant mass emitted within each class needs to be deter¬ 
mined. An analogous reasoning applies to area sources. 

A limitation of the method presented here is its use of a 
simple Gaussian plume model to calculate pollutant con¬ 
centrations within the short range (up to 100 km) around 
the emission source. The model does not consider local cli¬ 
matic effects, as they may occur in cities, for example. Fur¬ 
thermore, the model assumes a flat terrain around the source, 
i.e. the effect of topographical features such as valleys or 
hills, or of street canyons, is not taken into account. How¬ 
ever, the model is only used to calculate the area-integral in 
equation (3), rather than pollutant concentrations at indi¬ 
vidual points within that area. Errors in the calculated con¬ 
centrations at individual points are therefore cancelled out 
to some extent, or they are insignificant, if they only con¬ 
cern a small fraction of the entire area. For the latter reason, 
the effect of the street canyon in which the emission takes 
place turns out to be negligible, for example. While the street 
canyon increases the pollutant concentration by about 2-3 
times within a small area of some 10 meters in linear exten¬ 
sion, this increased contribution to <EE 1 near > is negligible 
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compared to the contributions from the much larger area 
within 100 km from the source (Nigge 2000) 4 . The extent 
to which <EE| ncar > is sensitive to built-up structures across 
the entire urban area (i.e. beyond the street canyon in which 
the emission takes place), to topographical features and to 
local climates, is a topic of further research. 

Given the continuous nature of the spatial variation of the 
population densities and meteorological variables, the in¬ 
troduction of discrete spatial classes of emission sites is in¬ 
evitably associated with a remaining variability of the im¬ 
pacts within the classes. A detailed statistical estimate showed 
this intra-class variability to be on the order of 40% of the 
class average impact. The averaging of the short-range ex¬ 
posure efficiencies over settlement structure classes and of 
the long-range exposure efficiencies over countries are con¬ 
sistent in their level of spatial detail (Nigge 2000). 

2 Site-Dependent Exposure Efficiencies 

2.1 Germany 

For the generic spatial classes in Germany, average exposure 
efficiencies <EE;> were calculated for emissions of a variety of 
primary pollutants (PM 2.5, PM 10, SO,, NO x , benzene, for¬ 
maldehyde, acetaldehyde, benzo[a]pyrene, 1,3-butadiene) at 
different effective emission heights between 5 m and 200 m 5 . 
They can be directly combined with exposure-response slopes 
such as the ones provided by Spadaro and Rabl (1999) or 
Hofstetter (1998) according to equation (1). Exemplary re¬ 
sults are presented in the following (for a complete list see 
Nigge 2000). An inhalation volume of 20 m 3 / (person x day) 
was used (Crettaz 2000). 

2.1.1 Traffic Emissions 

For emissions from ground-based transportation (cars, 
trucks, locomotives), an effective emission height of 5 m 
was assumed to account for a possible thermal plume rise 
or traffic induced turbulence. Given the fixed emission height, 
the focus in this section is on differences in the exposure 
efficiencies due to the influence of population density. 

The exposure efficiencies from traffic emissions of primary 
pollutants with short atmospheric residence times were found 
to vary by about a factor of 5 between large cities in 
agglomerations and rural districts in Germany if the annual 
mean wind speed is kept at a constant value (u = 3.5 m/s). 
This is shown in Fig. 4 for the example of acetaldehyde with 
an atmospheric residence time of 9 hours. Four out of the nine 
settlement structure classes are shown together with the coun¬ 
try average, which approximately represents the five remain¬ 
ing classes. Due to the short atmospheric residence time of 
acetaldehyde (9 hours), the exposure efficiencies are domi¬ 
nated by the short-range contribution. The spread by a factor 
of 5 in the exposure efficiencies is nevertheless smaller than 
the corresponding spread by about a factor of 20 in p 10 km (see 
Fig. 3). This is because EE near is dominated by contributions 


4 This refers to the effect of the street canyons on physical pollutant dispersion. 
Possible chemical effects of high background concentrations in street canyons 
on pollutant chemistry (e.g. of NO,) were not considered. 

5 A significant share of the total health impact of formaldehyde and benzo[a]pyrene 
is due to the uptake of food and drinking water (Hofstetter 1998), which is not 
considered here. 


from the first 10 km only for high values of p 10km (urban 
sites), while contributions from the range of 10 -100 km domi¬ 
nate for low values of p 10 km (average and rural sites). 

For substances with longer atmospheric residence times, the 
long-range contribution becomes more significant. This is 
shown for diesel particles (PM 2.5) with an atmospheric 
residence time of about 5 days in Fig. 5. Since the long- 
range contribution has the same value for all emission sites 
within one country, the relative spread between the settle¬ 
ment structure classes reduces to a factor of 2.2. At the same 
time, the absolute exposure efficiencies are increased for all 
settlement structure classes. 

The effect of a variation of the annual mean wind speed on 
the short-range exposure efficiencies is shown in Fig. 6 for 
the example of acetaldehyde. The country average exposure 
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Fig. 4: Spatial variation of exposure efficiencies EE from traffic emissions 
of acetaldehyde in Germany. EE D country average, p 10km population den¬ 
sity within the first 10 km, p 0 country average population density 
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Fig. 5: Spatial variation of exposure efficiencies EE from traffic emissions 
of diesel particles (PM 2.5) in Germany. EE 0 country average, p, okm popu¬ 
lation density within the first 10 km, p D country average population density 



u = 2.5 m/s u = 3.5 m/s u = 4.5 m/s 

Fig. 6: Influence of the annual mean wind speed u on the country average 
exposure efficiency EE Du from traffic emissions of acetaldehyde in Germany 


Int J LCA 6 (5) 2001 


261 



Generic Spatial Classes for Human Health Impacts, Part I 


LCA Methodology 


efficiencies vary by about a factor of 1.7 across the range of 
typical annual mean wind speeds in Germany 6 . In cases where 
the short-range contribution dominates the total exposure 
efficiency (i.e. for substances with short atmospheric resi¬ 
dence times, and for large cities in agglomerations even in 
the case of longer residence times), the annual mean wind 
speed therefore affects the total impacts to some extent. Its 
influence is nevertheless significantly smaller than that of 
the population density. Whether or not it is worthwhile to 
consider the influence of the annual mean wind speed can 
therefore be decided case by case, bearing in mind that this 
requires more detailed knowledge of the location of the 
emission source than a differentiation on the basis of the 
settlement structure classes alone. 


2.1.2 Influence of Emission Height 

The differentiation of exposure efficiencies between settlement 
structure classes decreases significantly with increasing emis¬ 
sion height. This is shown in Fig. 7 for emissions of particles 
(PM 10). The underlying effect is the decrease of the ground- 
level concentration within the first 10 km around the emis¬ 
sion site (Fig. 1). The urban settlement structure classes (top 3 
lines in Fig. 7) are most strongly affected by this decrease in 
concentration due to their high values of p 10 km . Exposure 
efficiencies for the other settlement structure classes (average 
and rural) are largely insensitive to emission height. The effect 
of the emission height h on the exposure efficiency, hence, 
depends on the settlement structure class. In the example of 
Fig. 7, a reduction of h by a factor of 20 (from 200 m to 10 m) 
increases EE by a factor of 2.2 for large cities in agglomerations, 
but only by a factor of 1.1 for low-density rural districts in 
rural regions. Vice versa, the influence of the settlement struc¬ 
ture class (approximately represented by p 10 km ) on EE depends 
on h, but is generally stronger: An increase of p 10 km by about 
a factor of 20 increases EE by a factor of 3.6 for h = 10m, but 
only by a factor of 1.7 for h = 200 m. 



d 
LU 
LU 
—. 

LU 
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Fig. 7: Influence of effective emission height and settlement structure class 
on the exposure efficiency EE from emissions of particles (PM 10) in Ger¬ 
many. Relative population density p, okm / p D indicated in brackets for each 
settlement structure class. EE 0 5m country average for h = 5 m 


6 The long-range contributions in Fig. 6 are averages across the zones of Ger¬ 
many to which the respective wind speeds apply. Their variation is geographic 
in nature rather than directly related to wind speed. 


Since five out of the nine settlement structure classes are 
always associated with an exposure efficiency very close to 
the country average, they can be merged into one generic 
spatial class. The spatial differentiation of the exposure 
efficiencies within Germany can therefore be captured by 
five generic spatial classes. Fig. 7 shows that the number of 
generic spatial classes can be further reduced for specific 
emission heights, e.g. to four classes for traffic emissions. 
For pragmatic purposes, it furthermore appears suitable to 
consider three typical effective emission heights: low (h = 
5 m, in particular traffic emissions), medium (h = 50 m) and 
high (h = 200 m). 


2.2 Europe 

The definition and operationalization of generic spatial 
classes to capture the spatial differentiation of health im¬ 
pacts of primary airborne pollutants with linear exposure 
response functions was demonstrated for the case of Ger¬ 
many. The extension of the generic spatial classes to other 
countries is a topic of further research. However, in order to 
obtain a first impression of the degree of spatial differentia¬ 
tion to be expected for Western and Central Europe, the 
average exposure efficiencies, EE country , for countries other 
than Germany were estimated. The Mediterranean Sea and 
the part of the Atlantic close to Europe were also consid¬ 
ered. The short-range contribution was approximately de¬ 
termined as 

^^'country, near ~ EE D ,near X (P counrr /Pn) x ( V u country ) (3) 

with 

Pcountry average population density for the respective coun¬ 
try (= 0 for emissions at sea) 

p D average population density for Germany (230 per¬ 
sons/km 2 , BBR1998a) 

u D annual mean wind speed for Germany (3.5 m/s, 
Gerth and Christoffer 1994) 

“country annual mean wind speed for the respective country, 
determined as a mean value of the measurement sta¬ 
tions for the country listed in the European Wind 
Atlas (Troen and Petersen 1989). 

The long-range contribution EE counffy far was calculated with 
the EcoSense software (IER 1998) as an average over sev¬ 
eral emission sites spread evenly over the area of the coun¬ 
try. Overall, 29 emission sites in Western and Central Eu¬ 
rope and 10 sites in the Mediterranean and the Atlantic were 
considered. For countries outside of the model area, which 
are relevant for the case study presented in Part II (Russia, 
OPEC countries), EE coumry>far was calculated as the product 
of a European fate factor determined as an average for the 
29 emission sites within the model area, an estimated effec¬ 
tive population density (9 persons / km 2 for Russia, 33 per¬ 
sons / km 2 as an average for the OPEC countries) and the 
inhalation volume V. 
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Fig. 8: Spatial variation of country average exposure efficiencies EE from 
traffic emissions of acetaldehyde across Europe (ATL Atlantic, MED Medi¬ 
terranean). EE d country average for Germany. Relative population density 
Pcountry 1 Pd indicated in brackets for each country 

Fig. 8 shows that the country average exposure efficiencies 
EE col ,ntry for traffic emissions of acetaldehyde vary by about 
a factor of 5.7/0.18 = 30. Most of the variation is due to 
Pcountry > while the influence of EE country hr and u count0 , is smaller. 
For diesel particles, the spread reduces to about a factor of 5 
due to their longer atmospheric residence time 7 . 

The country averages are likely to be similar to the values 
for emissions in rural areas of the respective countries (see 
Fig. 7). The high exposure efficiencies for traffic emissions 
in large cities in any country, on the other hand, can be ex¬ 
pected to be similar to those for Germany, i.e. on the order 
of about 3 EE d = 13.8 E-6 (Fig. 4). This is because they are 
largely determined by p, okm which is similar for large cities 
in all countries. The overall spread between the highest and 
the lowest class average exposure efficiencies across Europe 
can therefore be expected to be in the order of a factor of 
13.8/0.18 » 75 for primary pollutants with short atmospheric 
residence times such as acetaldehyde. An analogous estimate 
for diesel particles yields a lower spread by a factor of 8. 
Differentiations of this order can actually be expected to 
occur within individual countries that are scarcely popu¬ 
lated on average, but include large cities, such as within the 
Scandinavian countries. Generally speaking, the spread of 
the impact of traffic emissions between urban and rural ar¬ 
eas within one country can be expected to be the higher the 
lower the average population density of the country is. 

2.3 Comparison with other Studies 

Both the absolute levels of exposure efficiencies as well as their 
relative spatial differentiation are comparable with other studies 
from the literature. Potting (2000) found country average ex¬ 
posure efficiencies for emissions of hydrogen chloride at h = 
25 m across Europe to vary between 0.13 E-6 (Finland) and 
2.3 E-6 (Netherlands). Average exposure efficiencies for 
emissions of acetaldehyde (which has a similar atmospheric 
residence time) at h = 5 m in these countries were found to 
be about twice as high here (Finland: 0.27 E-6, Netherlands: 
4.6 E-6, Fig. 8), with a similar ratio between them. The dif¬ 


7 These factors also apply to other emission heights, since the country averages 
are not very sensitive to emission height (see Fig. 7). 


ferences in the absolute values are likely due to the different 
dispersion models used for the short range. The range of 
exposure efficiencies for PM 10 emissions in Germany be¬ 
tween 3.1 E-6 and 12.6 E-6 shown in Fig. 7 lies within a range 
across different emission heights and locations in France be¬ 
tween 1.7 E-6 and 15.7 E-6 for emissions of SO z (similar resi¬ 
dence time than PM10) as found by Spadaro and Rabl (1999). 
The larger range for France reflects its lower rural and higher 
maximum urban (Paris) population densities. European aver¬ 
age exposure efficiencies of 10 E-6 (Hofstetter 1998) for PM2.5 
and of 9.6 E-6 for fine particulates (Crettaz 2000) are compa¬ 
rable with the German average of 10.3 E-6 for the PM2.5 
emissions shown in Fig 5, but higher than the European aver¬ 
age of 6 E-6 determined in (Nigge 2000). The difference is 
probably due to differences in parameters or methods used 
for dispersion modeling. 

3 Conclusion 

This paper considered how the impact on human health of 
an emission of a primary airborne pollutant with a linear 
exposure-response function depends on the population den¬ 
sity around the emission site and on the emission height, 
and how this spatial differentiation can be considered within 
a Life Cycle Assessment. Concerning the dependence on the 
emission site, two spatial scales can be distinguished. On a 
larger scale, variations between countries or sub-continen¬ 
tal regions are due to their different average population den¬ 
sities. This sub-continental differentiation is insensitive to 
emission height. It can therefore be addressed by existing 
methods which are based on country average population 
densities and do not consider emission height. Across Eu¬ 
rope, the sub-continental differentiation amounts to about 
a factor of 30 for pollutants with short atmospheric resi¬ 
dence times of a few hours and about a factor of 5 for pol¬ 
lutants with long residence times of a few days. 

On a smaller spatial scale, there are also variations between 
urban and rural areas within one country. These deviations 
from the country average can be captured within LCAs by 
the method of generic spatial classes presented here. The 
strongest deviations from country average impacts apply to 
emissions at low heights in urban areas. Urban traffic emis¬ 
sions are a prominent example in this regard. It is for these 
situations that the presented method promises the highest 
benefit compared to methods based on country average popu¬ 
lation densities alone (Table 1). Traffic emissions are also 
the area where the application of generic spatial classes saves 


Table 1: Benefit of the proposed method over methods that consider na¬ 
tional averages of population density alone 8 



emission height 


low 

(<25 m) 

medium 
(25-100 m) 

high 

(>100 m) 

large cities in 
agglomerations 

++ 

+ | 

0 

'normal' cities b 

+ 

+ /0 

0 

average and rural 
districts 

0 

0 

0 


a scale: ++ / + / 0 (high / medium / low benefit) 

b cities in uibanized regions and highly densified districts in agglomerated regions 
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the most effort compared to a site-specific impact assess¬ 
ment, because the latter would have to consider a large 
number of individual sites along the path of the vehicle. 

For traffic emissions in Germany, the urban-rural differen¬ 
tiation was found to be about a factor of 5 for pollutants 
with short atmospheric residence times and about a factor 
of 2 for pollutants with longer residence times. It can be 
expected to be higher in countries with lower average popu¬ 
lation densities (i.e. most European countries). The need to 
consider the urban-rural differentiation in addition to the 
sub-continental spatial differentiation within LCAs there¬ 
fore increases with decreasing country average population 
density. Detailed calculations of site-dependent exposure 
efficiencies for generic spatial classes in countries other than 
Germany are a topic of further research. 

Three main limitations apply to the method of generic spa¬ 
tial classes presented here: It is based on the assumption of a 
linear exposure response function, which is valid for the 
major pollutants from energy generation and transportation. 
For other pollutants, further research needs to clarify to what 
extent spatial differentiation of impacts also arises from vari¬ 
ations in their background concentration. Further research 
is also required regarding the extent to which population 
exposures are sensitive to built-up structures, topographical 
features and local climates. The Gaussian plume model used 
here does not consider such effects. 

Inherent in the approach of generic spatial classes is a re¬ 
maining intra-class variability of impacts, which was found 
to be on the order of 40%. If more precise population expo¬ 
sures are required to reach an overall conclusion in any par¬ 
ticular application, a more detailed, site-specific assessment 
needs to be carried out. In such cases, the benefit of the 
method presented here lies in identifying the emission sources 
for which a site-specific assessment is worth the effort. 
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